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We consider a high density region of the hehum phase diagram, where the nuclei form a Bose- 
Einstein condensate rather than a classical plasma or a crystal. Helium in this phase may be present 
in helium-core white dwarfs. We show that in this regime there is a new gapless quasiparticle not 
previously noticed, arising when the constraints imposed by gauge symmetry are taken into account. 
The contribution of this quasiparticle to the specific heat of a white dwarf core turns out to be 
comparable in a range of temperatures to the contribution from the particle-hole excitations of the 
degenerate electrons. The specific heat in the condensed phase is two orders of magnitude smaller 
than in the uncondensed plasma phase, which is the ground state at higher temperatures, and four 
orders of magnitude smaller than the specific heat that an ion lattice would provide, if formed. Since 
the specific heat of the core is an important input for setting the rate of cooling of a white dwarf 
star, it may turn out that such a change in the thermal properties of the cores of helium white 
dwarfs has observable implications. 

I. INTRODUCTION 

The behavior of matter at high pressure is extremely rich, featuring the interplay of electromagnetic, statistical, 
quantum, and relativistic physics [1]. Aside from the intrinsic theoretical appeal of such a playground, the properties 
of matter in extreme conditions are important in many phenomenological applications. Perhaps foremost among these 
are applications in astrophysics and planetary science, as well as in terrestrial experiments on high-energy-density 
matter, for instance using inertial confinement or diamond anvil cells. Here we will discuss the behavior of helium 
under pressures high enough such that the average interparticle spacing / is smaller than the Bohr radius ag ~ 10^ fm. 
At such densities the helium gas becomes completely ionized, so that the system can be described as a plasma of 
helium nuclei and electrons interacting via electromagnetism. Such conditions are expected in helium-core white 
dwarfs. Helium white dwarfs (He WDs) are thought to be formed when a large part of a red giant star's envelope 
is removed prior to helium ignition, exposing the helium core[2-4]. A fair number of He WDs have recently been 
discovered[5, 6]. A particularly interesting aspect of these observations is the discovery that the sequence of the He 
WDs found in [6] comes to an early end compared to the sequence of carbon-oxygen white dwarfs. It may be that the 
explanation for this phenomenon involves developing a better understanding of the envelopes of the He WDs[6, 7], 
but it is also interesting to examine other alternatives, such as changes in the cooling rates of the WDs at lower 
temperatures due to phase transitions in the core. Hence our goal in this paper is to develop a better understanding 
of the ground state of the helium plasma in the cores of He WDs. 

Whether a plasma of nuclei and electrons is the ground state of a system depends on the density and temperature. 
There are two relevant possibilities for the behavior of high-density helium at lower temperatures if it cannot remain a 
standard plasma: crystallization or Bose-Einstein condensation of the helium nuclei. Which one is preferred depends 
on the density and temperature. For heavier white dwarfs, crystallization is the accepted model for low enough 
temperatures. The standard model for He WD evolution is based on assuming that the plasma of helium nuclei and 
electrons survives to low temperatures, see e.g. [8]. However, it was recently pointed out that in He WDs it is quite 
plausible that the conditions are such that the helium ions Bose-condense, with potentially observable consequences 
for the cooling rates of the stars[9-ll]. (For some early work on ion condensation in general see [12-18].) In this 
paper we will examine the consequences of helium ion condensation on the thermodynamic properties of the system, 
by developing an understanding of the spectrum of low-energy quasi-particle modes in the condensed phase. 

Let us now review the arguments that suggest that a helium plasma may Bose-condense at high enough densities. 
If the Coulomb energy Ec of the ions is much larger than their thermal energy Et and quantum effects can be 
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FIG. 1. (Color online.) Estimated crystallization (blue) and condensation (red) temperatures as a function of the WD core 
density, assuming the estimate for Tc in [11]. The crystallization temperature estimate in the plot assumes quantum effects 
are negligible, as discussed in Footnote 1. Bose-Einstein condensation of the nuclei is expected in the shaded area. 

neglected^, the ions should be expected to crystalhze. To get an estimate of the crystal-melting temperature, one can 
take Ec to be dominated by nearest-neighbor interactions, so that Ec — Z'^a/l, where Z is the ion charge, a ~ 1/137 
is the fine structure constant, and I is related to the density n through 3/(47r/'^) = n, while the thermal energy in the 
plasma phase is dominated by the contribution of the ions and can be approximated by that of a free bose gas Et ~ T. 
(We work with natural units where kg — c — h — \ except where stated otherwise.) The crystal-melting temperature 
Tmeit can then be estimated from the phenomenological relation Ec/Et ~ 180 [20-22] , giving Tmcit ~ (ao/07000K. 

On the other hand, if / becomes comparable to the thermal deBroglie wavelength Idb — \/27r/AfT, where M is 
the ion mass, one should expect the ions to Bose-condense. The critical temperature of such a nuclear Bose- 
Einstein condensate can be estimated by comparing this thermal wavelength to the inter-ion distance. This woud 
suggest a critical temperature of « 2tt/{MP) « 6.3/{MP). A free gas would condense at a smaller tempera- 
ture, Tc = [47rC(3/2)/3]-2/3 x 2Tr/{MP) « 1.27 /{Ml'^). Repulsive interactions, however, wih increase the critical 
temperature[23]. Some estimates[ll] suggest Tc « Att'^ /{3MP) « 13.2/(Af^^), a result qualitatively supported by 
more detailed calculations [24]. In any case, at a certain critical density, Tc will exceed Tmeit, and a Bose-condensed 
region of the phase diagram will open up. 

Numerically, assuming the estimate in [11] of Tc ~ 13.2/(M^^), the critical value of (ao/l) where the condensation 
region opens up is given by {ao/l) ~ 12.4, while using Tc ~ 6.3/{MP) gives the intersection value ao/l ~ 26. This 
translates to a critical mass density of around 2 x lO^'g/cm'^ to 2 x IQ^g/cm^, respectively. Meanwhile, helium white 
dwarf models show that their central mass density varies between 1 x 10^ g/cm^ to 1 x lO^g/cm'^, see e.g. [25, 26]. 
These estimates make it quite plausible that the conditions appropriate for Bose-condensation of helium nuclei may 
be present in the cores of Helium white dwarfs at some point during their thermal evolution. Despite (or perhaps 
because of) the uncertainties in the estimates, there is enough motivation to ask what the consequences of nuclear 
condensation may be for the bulk properties of a helium white dwarf, so that one can eventually decide whether the 
effects of nuclear condensation in He WDs may be detectable in current or future observations. The region of the 
phase diagram where Bose-Einstein condensation is likely to occur is shown in Fig. (1). 

At low temperatures, the thermodynamic properties of a system are determined by the spectrum of low-energy 
excitations. Any gapless or nearly-gapless quasiparticle modes make particularly important contributions to the 
thermodynamic properties, since the contributions of gapped quasiparticle modes become Boltzmann-suppressed once 
the temperature is low enough. Hence the first essential step in characterizing the physics of a quantum liquid is to 
understand the spectrum of quasiparticle excitations, and this is the focus of this paper. The most phenomenologically 
important result of our analysis is that the nuclear condensate has a previously undiscovered gapless quasiparticle 
mode. To find this mode it is crucial to pay careful attention to the constraints imposed by gauge invariance, and 
to characterize its properties correctly it is necessary to go beyond the Thomas-Fermi approximation in treating the 
electrons. Armed with the full quasiparticle spectrum, which turns out to be rather different in its details than 
appreciated in previous studies[9-ll, 27], we compute the specific heat of the helim-ion quantum liquid. We find that 
the specific heat of this quantum liquid is very small, far smaller than the specific heat of an ion crystal or a classical 
gas of ions. Fortunately, this means that our results do not qualitatively confiict with the promising phenomenological 



Quantum fluctuations may become important at high enough densities for very small temperatures, and can be expected to melt the 
lattice when the quantum zero-point energy of the ions Eq overwhelms the Coulomb energy keeping them localized. The density at 
which this happens is sensitive to the value that Ec / Eq must reach to trigger this quantum phase transition, and requires detailed 
calculation, see e.g. [19] in the case of hydrogen. Ref. [9] gives some plausible arguments that quantum effects melt the crystal at T = 
in He WDs, but fortunately for our purposes here it will not be important to take a position on the T = ground state in He WDs. 
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arguments regarding helium white dwarf coohng of Refs. [9-11]. Thus if the cores of hehum-rich white dwarfs undergo 
a phase transition to a phase with a nuclear condensate, one should expect a large drop in the specific heat of the 
cores, with a consequent increase in the rate of cooling of the stars. If this increase in the cooling rate is rapid enough, 
nuclear condensation may help explain the early termination of the He WD sequence observed by [6], but whether 
this is indeed the case is not yet clear as we discuss in the conclusions. 

The rest of this paper is organized as follows. Section II contains a description of the effective field theory describing 
the condensed phase, and presents a perturbative computation of the quasiparticle spectrum. Section III discusses the 
conditions under which corrections to the spectrum will be small, and is followed by a discussion of the implications 
of the quasiparticle spectrum for the specific heat of the condensed phase in Section IV. Finally, Section V contains 
the conclusions and a sketch of some directions for future work. 



II. EFFECTIVE THEORY AND QUASIPARTICLES 



Our main physics goal in this paper is to compute the specific heat of the nuclear condensate. The general principles 
of effective field theory and thermodynamics imply that the specific heat of a quantum liquid can be calculated from 
the dispersion relations of the low-energy quasiparticle excitations of the liquid. If the quantum liquid is the T ~ 
ground state, then the dispersion relations of quasi-particle excitations will be perturbatively close to the T = 
spectrum at low temperatures, so that in calculating the specific heat at small T one could use dispersion relations 
calculated at T = 0. In our case, however, the ground state of the system at T = may be a crystal, not a quantum 
liquid. In such a case, to compute the specific heat for Tmcit < T ^ Tc, one should use the quasi-particle spectrum 
computed at T = with the (in this case) artificial assumption that a condensate exists at T = 0. Temperature effects 
on the quasiparticle dispersion relations appear, together with other corrections, at loop level, and will be neglected 
here. Because the transition between the crystal and liquid phase (if there is one) must be first-order, there is no risk 
of finding unstable modes by perturbing around the condensed phase. Thus we will carry out our analysis at T = 
and assume that the helium nuclei are in a Bose-Einstcin condensate in order to compute the quasi-particle spectrum 
in the condensed phase and learn about its thermodynamic properties for T <^Tc. 

At the densities and temperatures of interest the typical interparticle distance / lies in the range Rn I ^ clq 
where i?jv ~ 2fm is the size of the a-particle, I is related to the density of ions n by ATrl^n/S = 1, and Qq = (q;to)~^ 
with m the mass of the electron is the Bohr radius. At these densities the helium nuclei can be treated as a point 
particles. An effective theory treating the ions, electrons and photons is described by the Lagrangian in the second 
quantized form 

C = ^^ [iD^ + A* + ^) ^ - \f^.uF^''' + X {iD^r + Me7" + m) X + • • • , (1) 

with Df^ip — ((9p — iZeA^)'i(j, Df^x = {dfj. + ieA^)x and Z = 2. The fields for the hefium nuclei and electrons are 
-0 and X respectively, and ji and /ig their chemical potentials. Omitted from eq. (1) are terms describing the strong 
force interaction between nuclei since they have a very short range, and Coulomb repulsion makes them irrelevant at 
the temperatures and densities we will be considering. 

To study the quasiparticle spectrum, we find it convenient to integrate out the electrons. The result is that the 
explicit fermion term above disappears and is replaced by a trace-log term coming from the fermion determinant: 

^electron = tr log (il?^7^ + ^i-el^ + m) 

- trlog {id^Y + A^e7" + H - ZenAo + \ j -^^A^^{-p)W^,{p)A%p) + 0{A% (2) 

where Zn is the electron density in an electrically neutral system. In the second line above we expanded the determi- 
nant in powers of the photon field. The first term above is of order 0{A^Y' and gives the free energy of a free electron 
gas. It makes an important contribution to the thermodynamics of the system and is associated with some gapless 
quasiparticles, the particle-hole excitations. However, this free contribution will not play a direct role in our discussion 
of the remaining excitations. The second term is the coupling between the electric potential and the electron charge 
that will cancel against the electrostatic energy of the ions. The third term above is of order 0{A^Y , and takes 
the form of a photon polarization tensor whose explicit form will be given below. Terms with higher powers of the 
photon field will not be needed in our discussion. The polarization tensor satisfies some relations as a consequence 
of gauge invariance. In particular, charge conservation implies that p^^^v{p) = 0- This relation, combined with the 
requirements of rotational invariance (but not boost invariance, which is broken by the electron chemical potential 
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term), implies that 11^^ is given by two scalar functions of po and p^, n(po,P^) and n''"(po,P^) as (see appendix) 



n 



p2 



(3) 



We will assume that the dynamics does not prevent the condensation of the bosonic ion field "01 in accord to 
the general discussion above, and parametrize as "0 = (^^ + /i)e"^. In order not to confuse excitations of physical 
degrees of freedom with the illusory excitations of gauge-dependent quantities we need to fix the gauge. We use the 
Fadeev-Popov method, choosing a variation of the i?^-gauge. The gauge-fixing term to be added to the Lagrangian 
is 



L 



gauge 



1 

2^ 



V.A doh-^——^ 

ev M 



(4) 



This gauge-fixing term cancels the quadratic terms mixing (j) with h or A. With our choice of field parametrization 
and gauge- fixing term, the Fadeev-Popov ghost fields are decoupled from h, (f>, and we can neglect them for our 
purpose of computing the spectrum. Up to quadratic order in the fields the Lagrangian becomes 



-Cquad = 2/it;ft. -I- Ze{v^ - n)Aa - l^i^^Y 



fih^ + 2ZevAoh 



2M' 



t,Zev 2M 



ZW^ 
2M ' 



2C 



(V-A)^ 



(5) 



where we are not showing the rest of the Lagrangian, £highor, which includes terms with more than two powers of the 
fields. Charge neutrality implies that {v + h)^ = n. We now choose unitary gauge by sending ^ — > oo. The advantage 
of unitary gauge is that it leaves only the "physical" fields in £quad- The main effect of the gauge fixing with this 
choice of ^ is to drop from the Lagrangian, as the term essentially gives it an infinite gap, causing it to 

decouple from the other fields. At the same time, going to unitary gauge eliminates terms involving doh from the 
Lagrangian. To proceed, we also separate into three parts: Aq, A^ and A^, with A(p) = A^(p) + A^(p), such 
that A-'-(p) • p = and A-'-(p) • A^(p) = 0. The quadratic part of the Lagrangian in unitary gauge is then given by 



Cquad = 2fj,vh + Ze{v^ - n)AQ + h 
I 



2M 



h + [-d^ + - -f- n^] A^ + Iao [-V^ + U] Aq 



8'^ 

•^0 



A^ + Ao 



dodi 



-1 



n 



Af+2ZevhA„, (6) 



with m\ = Z^e^w^/M = Z'^Airav^ /M . The chemical potential /i and the expectation value of (-0) = v must be 
chosen to guarantee charge neutrality which, at zero temperature, implies that — n. This requires fi — aX tree 
level. The physical interpretation of this choice is transparent: there is no Coulomb contribution to the energy from 
one photon exchange, since the ion and electron charge density cancel each other. The Coulomb energy arises from 
fluctuations of charge densities, described by two or more photon exchanges. These effects first arise only at one-loop 
order and are small. In [28] the one-loop value of the chemical potential was estimated to be of order /x ~ a/ 1, but 
we will not use this result here, only contributes to the propagator of the transverse photon A_l and it describes 
magnetic effects due to the electrons, such as Landau damping. Such effects do not play an important role here so we 
will simply neglect them. At quadratic order the two transverse modes decouple from /i,A^ and Aq, and have the 
dispersion relation (assuming is small) 



(7) 



The fields v4o, A^ and h mix so, in order to extract the spectrum of quasiparticles, we will integrate out two of them. 
The choice of which fields to eliminate is arbitrary. We integrate out A^ first. The propagator of A^ can be found 
through the easily verified formula 



{A5,j + Bd.djY 



1 

A 



(8) 



and find 



2M 



n^] A^ 



'^0 



n 



Aq + 2ZevhAo. (9) 
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Notice that if h and Aq had not coupled to A^, there would have been no time derivative in the h, Aq, sector and 
no quasiparticle could arise from this sector. However, A^ couples to h and Aq through the time-space component 
Iloj of the polarization tensor, which is required to be non-zero by gauge invariance if IIoo is non-zero. We can now 
integrate out either Aq or h. Eliminating we find 



Co 



2M 



2M- 



na^/v2 



n 



1 



A^ [- 



(10) 



The propagator Gh{p) for h is 



GhiPo,p) 



Pi 



PL 

2M 



p='+n(p,po) 



(11) 



The dispersion relations for the modes excited by h are given by the locations of the poles of Gh{po,p)- The fact that 
the residues at the poles have non-trivial momentum dependence has important consequences and will be discussed 
in Section III. It turns out that G/j has two poles describing (nearly) stable quasiparticles. Setting /i to its tree-level 
value n = 0, the positions of these poles are implicitly defined by 



Pq = 



9 2 



V2A/; p2+n(po,p)' 



(12) 



If the electrons were not dynamical, we would have 11 = and the mode would become gapped with po{p = 0) = ttia- 
This is a well known result [29-31] of the jellium model, which describes charged scalars interacting through an 
unscreened Coulomb potential in the presence of a static charge background to render the system electrically neutral. 

The combination of the generation of a gap for the transverse photons along with a massive scalar mode in the 
model when electrons are non-dynamical can be viewed is an example of the Anderson-Higgs mechanism. The electron 
sea, however, does contribute to the polarization tensor and alters the jellium model result. Already at the smallest 
density where the condensate is likely to appear gq/I « 12, the Fermi energy is about 0.5 MeV ^ 6 x 10^ K but the 
temperature is below 80 eV « 10^ K. Thus, the relevant temperatures are far below the Fermi energy, and we can use 
the zero temperature polarization tensor. In general, n(po, p) has a complicated non-analytic dependence on pq and 
p. The behavior of TI{pq, p) as po and p become small, the relevant regime for the low lying excitations we are after, 
depends on how the limit is approached. The asymptotic forms of 11 that we will need can be found in e.g. [32-34] 



lim Ren(p,po = 0) = lim '^^^<^^P' 
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1-a; 



1 -p/2kp 
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(13) 



where fie = \/k'^ + m? is the electron chemical potential, vp = kp/fie is the Fermi velocity. The non-relativistic 
result valid at smaller densities is obtained from the above by setting — > to and vp — )■ kp/m. 

By definition gapped quasiparticles have a finite energy po as their momentum approaches zero. In this regime 11 
is given by 



n « - 



AnaZn p^ 



fJ-e 



Pi 



(14) 



In the p limit, eq. (12) gives a pole at 



pI 



e'^Zn 



= AirZan 



M 



1 

Me 
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^red 



(15) 
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with rrired = fJ-eM/^Zfie + M) and Zn = rig is the electron density. This mode is gapped, justifying the assumption 
that po remains finite as p — 0. Its frequency is changed from the electron plasma frequency ujp = ATrarie/rn by (i) 
relativistic effects changing m into /Xg and (ii) the change of the electron mass to the reduced mass, since the electrons 
are oscillating against the nuclei and not against an static positive background. Numerically we find 



Po 



500 eV 1 



^-2 



5.8 X 10" K 1 



-1/4 



oo 
I 

-1/4 



/ ■ 



(16) 



We can also look for gapless modes by assuming that po — xvpp {x ^constant) as p — > 0. Wc have 

Ren^m^ (^l_|log|l±||) 
Imll — —m1xQ{l — x) 
with = Aa^ekp/"^- As p the p** term in eq.(12) can be neglected, and we find 

Po _ 2 2 _ ™A 
~t: — X Up — 



(17) 



(18) 



(19) 



The left-hand side is positive for x < 0.8355 and reaches a maximum of ss 0.211 at a; = 0.623. The association of 
the solution of eq. (19) with a quasiparticle implicitly assumes that the imaginary part of the quasi-particle's energy 
(that is, the attenuation constant) is small, which as we will see would not be the case for large x. Numerically, the 



or 





l + x 




(l-|log 
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parameter 



is given by 



rriA 
vprUs 



0.8 ^2/3(1 



0.0012(1 




(20) 



and, to a very good approximation x ~ niA/ {vp'^s) and is very small. So eq. (19) can be solved. The solution of 
eq. (12) with IT given by eq. (17) can be found using an expansion in powers of the small parameter a;, and is given by 



rriA 
Po = P 



(21) 



The imaginary part is suppressed compared to the real part by a factor of ■Kx/{Avf) <C 1, and this gapless mode 
is almost stable. If we had used the Thomas-Fermi approximation for the screening of the electrons, 11 = m^, the 
imaginary part of the dispersion relation would not have been visible. The imaginary part encodes the effects of 
inelastic scattering of this mode with particle-hole excitations of the electrons. The size of the imaginary part is 
actually enhanced by the high density parameter a^/l, but this enhancement is overwhelmed by the suppression 
coming from the mismatch between the rest masses of the ion and the electron. It seems reasonable to think that the 
physical reason for the near-stability of the gapless mode is that collisions between it and the particle-hole excitations 
are mostly elastic thanks to the large mismatch between the masses of the helium nuclei and electrons. 

This gapless mode was not noticed in previous analyses of the system. The inclusion of the off-diagonal elements 
proportional to Hon in the polarization tensor, demanded by gauge invariance, are crucial for their existence. One 
interpretation of this mode may be as some kind of electron zero-sound-like excitation mixed with nuclear condensate 
excitations, but its velocity and attenuation rate behave differently than those of zero sound in Fermi liquids. As we 
have already noted this mode would not survive if the electrons were not dynamical. It would also not survive if the 
nuclei were not dynamical and the condensate were to vanish, since as is well-known an electron gas in a neutralizing 
positive non-dynamical background does not have a zero-sound-like mode - instead there are only the (gapped) plasma 
oscillations. In section III, where we discuss the possible corrections to the dispersion relations, we will have more to 
say on the identity of this gapless mode. 



7 



The linear dispersion relation is valid for all relevant momenta at the temperatures we are interested in. In fact, 
the (real part of the) polarization 11, including the next order in the momentum expansion, has the form 



n(po 



xvfP,P) 



Pi 



9{x) 



(22) 



where g{x) is a function of order one, assuming II{pq = xvpp, p) is analytic in around p^ = 0. Then the form of 
n we show above follows essentially from dimensional analysis, since the only relevant momentum scale in the loops 
defining 11 is kp. The solution to eq. (12) is 



771?. 



m1g(x) 



1 



(23) 



The p'^ term becomes as important as the p^ term for p = p : 



The energy of the gapless mode at this 
momentum is po = toa/'tIsP « 77ia- The scale 777 A^ however, is significantly higher than the temperatures where the 
condensate can exist. Numerically it is given by 



TUA 



= \/iZa^7 



1.3 X 10'' K 



(?) 



3/2 



'ao\3/2 



(24) 



which translates to a temperature of 2.6 x IQ^K for an He WD with a central mass density of 5 x lO^g/cm'^, or 
Qq/I ^ 35, which we will take as our reference density parameter for white dwarfs. 



III. HIGHER ORDER CORRECTIONS 



Up to now our analysis was restricted to leading/tree level order. We will now discuss higher order corrections in 
order to verify that our leading results are indeed reliable. This is specially important for the gapless mode as no 
apparent symmetry exists that would protect this mode from picking up a gap when corrections are calculated. As 
was the case in the context of the deuteron condensate discussed in [28] , we expect that the perturbative corrections 
to our results above will be in powers of ao/l. However, the existence of a gap for the gapless mode we found at tree 
level, even if small, could qualitatively change our conclusions. 

One way of organizing the perturbative expansion for our model is to integrate out all but the h and A-^ fields, 
as it was done in eq. (10), but keep terms with higher powers of h. The residue of the propagator for the h field, 
Gh{p) in eq. (11) has a non-trivial momentum structure. For this reason it is more convenient to work with the 
canonically-normalized field H defined hy H = {—V^ /M)^^/^h whose propagator is 



Gh{po,p) 



1 



1 



pI 



Pi 

2M 



P "ij 



hn(p,po) 



Pi 



(25) 



hp-' 



where the approximation above is valid near the propagator's pole for small momenta. The crucial point is to observe 
that the self- interactions of H described by terms proportional to H^, H'^, ■ ■ ■ always include derivatives. As a result, 
the action for the H field has the symmetry H — > H + 7], with 77 a spacetime constant. This symmetry guarantees that 
a gap will not be generated at any order in perturbation theory (or even non-perturbatively) . In fact, this is the same 
symmetry that guarantees that Nambu-Goldstone bosons remain gapless to all orders in the loop expansion. For all 
practical matters, H seems to be a Nambu-Goldstone boson, despite the fact that the symmetry being spontaneously 
broken, the electromagnetic C/(l), is a local symmetry. Before commenting on this somewhat surprising result, let us 
verify that H is indeed derivatively coupled. 

We claim that a term involving a number 7i oi H fields contains at least n derivatives, that is, they are schematically 
of the form with m > n, where d stands for a time or space derivative. If this statement were false, there 

would be terms in the effective action of h of the form Q™'-"H^ with m > n, and we now argue that such terms are 
impossible. In fact, self-interaction terms for h can be separated in two groups. One consists of terms generated by 
integrating out the fields Aq and A^. An example is shown in Fig. (2). These diagrams, at small external momenta, 
generate terms proportional to a/m'ld'^H^. Notice that, if the ^0 photons were not screened, diagrams like those in 
the left panel of Fig. (2) would generate terms exactly of the form we claim is impossible. The screening of Aq photons 



FIG. 2. (Color online.) Some of the graphs contributing to the self-interactions of the h mode. Dotted lines represent the h 
propagators, wavy lines are Ao photons, while solid lines are electrons. 



is essential to this result. Another class of contributions is encoded in diagrams arising from the terms of order 
and higher in eq. (2), which describe interactions mediated by the polarization of the electron sea. An example is 
shown in the right panel of Fig. (2). Again, if the photons were massless, these diagrams would generate terms of 
the kind we claim is forbidden, but the screening of the Coulomb force eliminates this possibility. In addition, it is 
conceivable that the non-localities of the fermion loop might generate negative powers of the external momenta. It is 
easy to see that this is not the case. In fact, a fermion loop with n photon legs attached to them is proportional to " 

) ^^Lll^fd^k \, ^—j Zn ^m3/V^/2-"+i (26) 

fco + Me-H+iOfco; d^^^^ J ko + ^ie - + lOko d^Ji^\^^^ 

density 

(for simplicity we assumed non-relativistic electrons in showing this expression). The essential point is that the loops 
are infrared finite and the diagram does not diverge as the external momenta vanishes'^. 

The conclusion we derive from this analysis is that, for all practical effects, the gapless mode we found is indistin- 
guishable from a Nambu-Goldstone mode. This may seem at odds with the fact that the only symmetry which is 
spontaneously broken, the U{1) gauge invariance, is actually a local symmetry. In relativistically invariant theories 
one would expect the Anderson-Higgs mechanism to play its usual role, and no physical gapless particle to arise if the 
symmetry being broken is local, but the situation is more subtle in non-relativistic theories. In the non-relativistic 
case, it turns out that much rests on whether the electric forces mediated by the photons are screened. Recall that 
the Anderson-Higgs mechanism avoids the conclusion of the Goldstone theorem due to infrared divergences resulting 
from long-range forces mediated by unscreened photons. In our model, the screening due to the electrons destroys 
the long range forces, and the Goldstone theorem remains in effect. This phenomenon is actually well known, for a 
general discussion see ref. [35-37]. In a model where there is screening but no dynamical electrons, there would be a 
gapless mode, stable at tree level, as can be easily seen from eq. (12) with 11 = m^. In our model the gapless mode 
can transfer energy to the electrons and decay already at tree level. Due to the difference in masses between ions and 
electrons this transfer is inefficient and the mode is almost stable. 



IV. SPECIFIC HEAT 



As pointed out in Refs. [9-11] the main observable impact of having a nuclear condensate instead of an ion lattice 
or an ion gas inside white dwarfs is a change in the cooling curves caused by a change in the specific heat. The 
specific heat from the condensed phase at low temperatures comes both from the electrons and the new gapless mode. 
The contribution of the plasmon mode and the massive transverse photons is Boltzmann-suppressed. The electron 
contribution to the specific heat per ion originates from the first term in eq. (2) and is well known: 



(27) 




^ Incidentally, this result in the n = 2 case reproduces the value of ^ mkp. 

^ The value of these loops depend, at small momenta, on the ratio po/p- But for all values of po/p the zero-momentum limit is finite. 
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FIG. 3. (Color online.) Specific heat per ion from electrons (red) and nuclei (blue), with the temperature is held fixed at 
T = 5 X 10^ K. The band indicates the regions with densities relevant for He WDs, which range between 3 x 10* gr/cm^ and 
10® g/cm? . The corresponding values range from 5.2 x IQ^K to ~ 1.2 x IQ^K. Also shown is the value for free ions (3/2 in 
yellow) and lattice phonons that would exist if the nuclei crystallized (dashed lines). 



The gapless zero-sound-likc mode contribution can be easily calculated through 

-T^ (28) 



I de 1 d f d^p chP 1 Stt^ ^3 



" ndT ndT J (27r)3 e^P/^ - 1 n 15c|^ ' 
for a mode with dispersion relation po = chP- Numerically we have 

16^5/2^3/2 / ; \ 9/2 / J. 



Cy 



15^/3 \m J \ao 



m 



2 / ; \ 9/2 



. 3.3 X 10^1 + ^)-^-J ^—J . (29) 

The contribution of the gapless NG-likc mode to the specific heat obviously has a very different temperature 
and density dependence than the contribution of the electrons, and which one is dominant depends on the specific 
parameters of interest. In Fig. (3) we show the values of the electron and ion to the specific heats for the fixed 
temperature T = 5 x 10^ K, at densities relevant for white dwarfs. For this particular choice of parameters, the 
electron contribution happens to be dominant, but at higher temperatures the situation can be different. If instead 
of an nuclear condensate we had an ion lattice, as might be the case at low enough temperatures, the specific heat 
would be dominated by the lattice phonon contribution. This contribution has been analyzed before [38] and is given 

by 



lattice 



167r4 / T ^ ^ 



leTT* 1 /M\^/^ f iV^ f T 



1.5 X 10« f^V^V-— - ) , (30) 



aaj \10^ K 

where the Debye temperature 0d is the same as what we called niA, and the formula given above is valid for T ^ Ojj. 
This contribution is enhanced, as compared to the condensed phase, by a power of The specific heats from 

the different sources are plotted in Fig. (3) . The picture for the specific heat of the nuclear condensate we have found 
broadly agrees with the results of [9-11]. Taking into account the contributions of all relevant quasiparticles, we see 
that if nuclear condensation takes place in the cores He WDs, their specific heat of the cores will be much lower than 
if the helium were in a crystalline or plasma phase. 



V. CONCLUSIONS 



The quasiparticle spectrum of the nuclear condensate turns out to be unexpectedly rich. In addition to transverse 
photons picking up a magnetic mass through the Higgs mechanism and the existence of a gapped plasmon quasiparticle. 
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we also found that the system supports an almost stable gapless mode. Somewhat surprisingly, despite the fact that 
the symmetry being broken is a local one, this physical gapless mode can be identified as a Nambu-Goldstone mode, 
arising because in our model the condensing field is non-relativistic and because the electrons screen the long-range 
Coulomb forces. To find this mode it turned out to be crucial to take into account the constraints of gauge invariance 
on the electron contribution to the photon polarization tensor. This gapless excitation has a number of interesting 
properties from the perspective of field theory, and their further exploration is a promising direction for future work. 
For phenomenology it is important to explicitly compute the finite-T corrections to the spectrum, since the relevant 
temperatures for He WDs are not so widely separated from Tc, and finite-T corrections could be important for some 
quantities. An effect we did not consider in this paper is that at low enough temperatures the electrons will pair 
up and condense, and it would be nice to understand better when such phenomena become relevant. It would also 
be interesting to check whether there is an analogous mode in super-dense deuterium [28, 39], as seems likely, and to 
understand its phenomenological consequences. 

With the spectrum of quasiparticle excitations in hand, we briefly examined the specific heat of the nuclear con- 
densate, comparing our results to previous work in the literature. For the parameter values appropriate for describing 
helium white dwarfs, the contribution of the gapless (Nambu-Goldstone-like) mode to the specific heat can be larger 
or smaller than the contributions from the particle-hole excitations of the degenerate electron gas, but generically both 
contributions are very small compared to that of an uncondensed plasma. Our results suggest that the specific heat 
of the core of a white dwarf will drop by a factor of ~ 10~^ upon ion condensation, similarly to previous analyses [10]. 

It remains an interesting open question whether such a change in the thermodynamic properties of the core of a 
white dwarf would have observable consequences. This matter was previously discussed in Ref. [10], which argued that 
the consequences would be detectable based on a simple model, and in Ref. [26], which argued that observable effects 
were unlikely based on a somewhat more detailed model. The ultimate answer to this question depends on developing 
a detailed model of white dwarf cooling in the presence of ion condensation, taking into account all relevant features 
of the system, such as e.g. the phenomenological consequences of the suggestion of Ref. [24] that the transition to the 
condensed phase is first order rather than second order, or the implications of the complete quasi-particle spectrum 
on e.g. the rate of neutrino emission[40]. Aside from looking at the effects of condensation on the cooling rates of He 
WDs, it may also be interesting to consider whether the way electric and magnetic fields are screened in the nuclear 
condensate may have any observable consequences for He WDs[41, 42]. 
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VI. APPENDIX: POLARIZATION TENSOR IN THE ABSENCE OF BOOST INVARIANCE 

In the main text we used relations between the different components of the polarization tensor H^jy imposed by 
gauge symmetry even in the absence of boost invariance. In particular we used the fact that, in the absence of boost 
invariance, H^i, depends on two arbitrary functions of po and — PiP^ , i = 1,2, 3. In order to see this, consider the 
most general symmetric tensor, function of p^ and (where — (1, 0, 0, 0) points in the time direction): 
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n, 



A{po,p'^)p,,Pu + B{po,p^)nf,n^ + C{po,p'^){pf,n^ + p^Uf,) + D{po,p'^)gi 



(31) 



with pq = p^n'^ . The Ward identity implies in 



= p^H 



A{P(i^P^)P^Pu + B{pQ,p'^)pan^ + C{po,p'^){p'^n^ + P^Pq) + D{pq,p'^)p, 



(32) 



This means that we can eliminate two of the four arbitrary functions using the relations 



^(po,p')/ + C{po,p^)pQ + D{po,p^) = 0, 
B{po,p^)po + C{po,p^)p' = 0. 



(33) 
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